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ABSTRACT 



We monitored the chromospheric activity in the Ca II H & K lines of 13 
solar- type stars (including the Sun); 8 of them over three years at the CFHT 
and 5 in a single run at the VLT. Ten of the 13 targets have close planetary 
companions. All of the stars observed at the CFHT show long-term (months to 
years) changes in H & K intensity levels. Four stars display short-term (days) 
cyclical activity. For two, HD 73256 and Ceti, the activity is likely associated 
with an active region rotating with the star, however, the flaring in excess of 
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the rotational modulation may be associated with a hot jupiter. A planetary 
companion remains a possibility for Ceti. For the other two, HD 179949 and 
V And, the cyclic variation is synchronized to the hot Jupiter's orbit. For both 
stars this synchronicity with the orbit is clearly seen in two out of three epochs. 
The effect is only marginal in the third epoch at which the seasonal level of 
chromospheric activity had changed for both stars. Short-term chromospheric 
activity appears weakly dependent on the mean K-line reversal intensities for the 
sample of 13 stars. Also, a suggestive correlation exists between this activity 
and the MpSini of the star's hot jupiter. Because of their small separation (< 
0.1 AU), many of the hot jupiters lie within the Alfven radius of their host stars 
which allows a direct magnetic interaction with the stellar surface. Wc discuss 
the conditions under which a planet's magnetic field might induce activity on the 
stellar surface and why no such effect was seen for the prime candidate, r Boo. 
This work opens up the possibility of characterizing planet-star interactions, with 
implications for extrasolar planet magnetic fields and the energy contribution to 
stellar atmospheres. 

Subject headings: stars: late-type, activity, chromospheres, planetary systems, 
radiation mechanism: non-thermal, stars: individual: r Boo, HD 179949, HD 
209458, 51 Peg, v And, HD 46375, HD 73256, HD 75289, HD 76700, HD 83443, 
Get, T Get 



1. Introduction 

In the past decade, ^ 130 extrasolar planets have been discovered as companions to 
late-type dwarfs mostly from Precise Radial Velocities (PRV; e.g. Walker et al. 1995) from 
which the Keplerian orbital parameters and minimum planetary masses can be derived. 
From transiting planets, of which five are known, true planetary masses, radii, densities 
and albedo-dependent estimates of surface temperature are available. To date, a handful of 
experiments have lead to further information about these planets, including Lyman a (Vidal- 
Madjar et al. 2003) and sodium detections (Gharbonneau et al. 2002) in the atmosphere of 



^ Based on observations collected at the Canada- Prance-Hawaii Telescope operated by the National Re- 
search Council of Canada, the Centre National de la Recherche Scientifique of France, and the University of 
Hawaii as well as data from the European Southern Observatory's Very Large Telescope, Chile (programme 
ESO No. 73.C-0694). 



-3- 



the transiting planet around HD 209458. A current review of follow-up techniques to probe 
planetary characteristics has been presented by Charbonneau (2003). 

The presence of a giant planet likely influences its parent star beyond the dynamical 
perturbations measured by the PRV method. Robinson & Bopp (1987) caught the signa- 
tures of "superflares" with energies of ~ 10^^ ergs (Schaefer et al. 2000) on nine solar analogs 
which have no otherwise unusual properties such as very rapid rotation or very high chromo- 
spheric activity. Rubenstein & Schaefer (2000) suggested that these anomalous superflares 
were stimulated by an unseen close-in extrasolar giant planet (CEGP, also known as a 'hot 
jupiter'). They explored the possibility of magnetic reconnection events occurring between 
both the star's and the planet's entangled magnetic fields. 

Cuntz et al. (2000) suggested a more consistent observable interaction between a parent 
star and its hot jupiter in the form of enhanced stellar activity of the star's outer atmosphere. 
This interaction can take the form of tidal and/or magnetic heating of the plasma. Both the 
star and its planet experience strong tidal forces as well as repeated magnetic reconnection 
with a large planetary magnetosphere. If such planet-induced heating of the star is confined 
to a narrow range in stellar longitude, the heated regions likely track the planet as it orbits. 
This implies that the period of any observed activity would be correlated with the planet's 
orbital period Porb such that tidally induced activity has a period of ~ Porbf^ and magnetic 
activity, a period of Porb- In the simplest configuration of magnetic interaction, the expected 
enhancement would occur near the sub-planetary point, which defines orbital phase = 
at the time when the planet is in front of the star relative to the line-of- sight. 

A partial analogy can be made to the Jupiter-Io system (Zarka et al. 2001) where the 
volcanically active moon of Jupiter, orbiting at a distance of 5.9 Rj, constantly couples with 
Jupiter's magnetosphere, leaving two footprints at high positive and negative latitudes on 
the planet's surface. Plasma flows along the magnetic field lines, making up the lo Flux 
Tube, as these footprints follow lo in its orbit. Even though the analogy is limited, a similar 
phenomenon may occur between CEGPs and their stars such that coupling between the 
magnetic fields of the planet and the star may cause footprints or "hotspots" on the star's 
surface which follow the planet and have a period close to the planet's orbital period Similarly, 
auroral emission may be stimulated on the planet's atmosphere (Zarka et al. 2001) but no 
searches for planetary radio emissions have yet been successful (e.g. Bastian et al. 2000, 
Farrell et al. 2003, Ryabov et al. 2003). 

More recently, Santos et al. (2003) observed cychc photometric variations with a period 
very similar to that of the radial velocity (RV) curve for the K dwarf HD 192263. The 
stability of the 24.4-day periodic radial velocity through almost 4 years of data (Santos et 
al. 2000, Henry et al. 2002) rules out the interpretation that stellar activity alone is the cause 
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of the RV curve and supports the existence of a planetary-sized companion around the star. 
However, they question what might cause a quasi-stable photometric period that coincides 
with the planetary orbit. They offered planet-induced magnetic activity off-set by 90° from 
the sub-planetary point as an explanation. 

Even though the interpretations of Rubenstein & Schaefer (2000) and Santos et al. 
(2003) arc uncertain, there exists ample observational evidence of such tidal and magnetic 
interactions in the exaggerated case of the RS Canum Vcnaticorum (RS CVn) stars, which 
are tightly-orbiting binary systems consisting of two chromospherically active late-type stars. 
For example, Catalano ct al. (1996) found as many starspots and plages within 45° of the 
sub-binary point^ as on the rest of the stellar surface for several RS CVn systems. A And, a 
relatively long-period system {Porb — 20.1 d; Walker 1944), shows modulation of the Mg II 
UV chromospheric emission lines with a period of 10 days, half the orbital period. Glebocki 
et al. (1986) interpreted this as a tidal heating of the primary by its companion, possibly a 
brown dwarf (Donati et al. 1995). And, in our own Ca II H & K observations of ER Vul, 
an RS CVn system with two G V dwarfs and Porb — 0.69 d, we see clear enhancements near 
the sub- binary longitudes (Shkolnik et al. 2004). 

With these scenarios in mind, we searched for periodic chromospheric heating by mon- 
itoring the Ca II H & K emission in stars with giant planets within a few stellar radii. 
We chose to study the tightest observable systems since the tidal and magnetic interactions 
depend on the distance from the planet to the star as and respectively. Of the 

known extrasolar planets, about 20% have semi-major axes of less than 0.1 AU and masses 
comparable to Jupiter's (Schneider 2004). It is expected that these planets have magnetic 
fields also similar to Jupiter's (4.3 G).^ It is also reasonable to assume that any magnetic 
interaction would be greatest in the outermost layers of the star, namely the chromosphere, 
transition region and the corona due to their proximity to the planet, low density, and non- 
radiative heat sources. The broad, deep photospheric absorption of the Ca II H & K lines 
allows the chromospheric emission to be seen at higher contrast. Because of this and the 
accessibility from the ground, the H & K reversals are an optimal choice with which to 
monitor chromospheric heating of these sun-like stars. 

Our program stars have orbital periods between 2.5 and 4.6 days, eccentricities ~ 
and semi-major axes < 0.06 AU. These systems offer the best chance of observing upper 
atmospheric heating. The first five systems we observed were r Boo, HD 179949, HD 209458, 



^The 'sub-binary' point refers to the longitude on the star that faces the binary system's center-of-mass. 

^However, if tidally locked, the planets' rotation rates may be much lower than Jupiter's possibly causing 
their magnetic fields to be substantially smaller (e.g. Sanchez-Lavega 2004). 
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51 Peg and v And from the Canada- France-Hawaii Telescope (CFHT). The first results from 
2001 and 2002 observations, including the first evidence of planet-induced magnetic heating 
of HD 179949, were published in Shkolnik et al. (2003). We later extended the experiment 
at the Very Large Telescope (VLT) to include five southern targets, HD 46357, HD 73256, 
HD 75289, HD 76700, and HD 83443. The system parameters for the program stars are 
listed in Table 1 along with our two standards, r Ceti and the Sun. 

In this paper, we compile our 2003 CFHT data with those of previous years and include 
the recent VLT observations. A broader understanding of stellar activity, its cycles, and 
planet- induced chromospheric heating emerges. We also observed Ceti, a young (650 — 
750 Myr old), chromospherically active, solar analog. It was one of the nine stars reported to 
have anomalous superflare activity by Robinson & Bopp (1987) possibly caused by magnetic 
interactions with a close-in giant planet (Rubenstein & Schaefer 2000). Of their sample, only 
k} Ceti has been looked at by the PRV planet searches of Walker et al. (1995), Cumming et 
al. (1999) and Halbwachs et al. (2003); none of which has detected a planet. Interestingly 
though, Walker et al. (1995) observed a rapid RV change of 80 m s~^ in 1988 which could 
have been caused by a planet in highly elliptical orbit. However, it coincided with an equally 
sharp increase in the Ca II chromospheric activity indicator at 8662 A implying the RV jump 
was likely due to changes intrinsic to the star. We may have observed this same phenomenon 
in our Ca II data. 

The details of our CFHT and VLT observations arc outlined in Section 2. We briefly 
discuss the precise differential radial velocities which yielded updated cphcmcridcs for the 
planetary systems such that orbital phases could be determined to better than 0.02. In 
Section 3 we discuss our analysis and results of the Ca II H & K measurements including 
long-term, short-term and rotational modulation. A theoretical discussion of the physical 
requirements for magnetic interactions between stars and their hot jupiters is presented in 
Section 4 with future experiments suggested in Section 5. 
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Table 1. Stellar and Orbital Parameters 



Star 


Spectral 


vs'mi 


Prot 


a b 




c • • b 

Semi-major 




(K') 


(MAD K)** 




Type 


(km s ) 


(days) 


(days) 


{M.j) 


axis (AU) 


(A) 


(A) 


(A) 


T Boo 


F7 IV 


14.8±0.3 


3.2* 


3.313 


3.87 


0.046 


0.323 


0.177 


0020 


HD 179949 


F8 V 


6.3±0.9 


<9S 


3.092 


0.98 


0.045 


0.399 


0.202 


0.0051 


HD 209458 


GOV 


4.2±0.5 


16^ 


3.525 


0.69' 


0.045 


0.192 


0.077 


0.0009 


51 Peg 


G2 IV 


2.4±0.3 


21.9* 


4.231 


0.47 


0.05 


0.178 


0.071 


0.0008 


V And 


F7 V 


9.0±0.4 


14f 


4.618 


0.7lj 


0.059 


0.254 


0.091 


0.0016 


HD 46375 


Kl IV 


<2 


43"^ 


3.024 


0.25 


0.041 


0.339 


0.233 


0.0011' 


HD 73256 


G8/K0 V 


3.22±0.32 


14™ 


2.549 


1.85 


0.037 


0.995 


0.899 


0.0051" 


HD 75289 


GO V 


4.37 


15.95™ 


3.510 


0.42 


0.046 


0.177 


0.062 


0.0003 


HD 76700 


G6V 


<2 




3.971 


0.20 


0.049 


0.212 


0.095 


0.0004 


HD 83443 


KOV 


1.4 


35"^ 


2.986 


0.38 


0.039 


0.314 


0.211 


0.0010 


Get 


G5 V 


4.64±0.11 


9.3° 








0.906 


0.815 


0.0059" 


r Get 


G8 V 


1 


33 








0.201 


0.102 


0.0003 


Sun 


G2 V 


1.73±0.3 


27 








0.298 


0.137 


0.0005 



*The first five periods are from this work. (See Section 2.1.1.) See Note b for the last five. 

''PubUshed orbital solutions: r Boo & v And - Butler et al. 1997; HD 179949 - Tinney et al. 2001: 51 Peg - 
Marcy et al. 1996; HD 209458 - Gharbomicau et al. 1999); HD 46375 - Marcy et al. 2000; HD 73256 - Udry et 
al. 2003; HD 75289 - Udry et al. 2000; HD 76700 - Tinney et al. 2002; HD 83443 - Mayor et al. 2004. 

'^Total integrated intensity bounded by the Kl features of the mean normalized Ga II K core. These values are 
relative to the normalization points near 3930 and 3937 A which are at approximately 1/3 of the continuum at 3950 

A. 

''We subtracted the photospheric emission from (K) in order to measure the mean integrated chromospheric 
emission (K') using data from Wright et al. (2004). (See text for more details.) 

''Average integrated 'intensity' of the mean absolute deviation (MAD) of the K residuals, per observing run 

*Henry et al. 2000 

^Prot is calculated from the vsini and a stellar radius of 1.24 Rq (Tinney et al. 2001). 

■^Mazch ct al. 2000 

'Transiting system; i = 86.1° ± 0.1° (Mazch ct al. 2000) 
^Glosest of three known planets in the system 

''Derived from empirical fits (Noyes et al. 1984) of the Ga II iJ^^ index (Wright et al. 2004) 

'(MADK) was corrected for a « 30% variable contribution by telluric Ga II emission since it was the first star 
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observed after sunset. 
™Udry et al. 2003 

"Value corrected to remove modulation due to rotation. For HD 73256, the non-corrected value is 0.0155 and for 
Ceti, 0.0182. 

°Rucinski et al. 2004 
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2. The Spectra 

2.1. CFHT Observations 

The observations were made with the 3.6-m Canada-Prance-Hawaii Telescope (CFHT) 
on 3.5, 4, 2, and 5 nights in 2001 August, 2002 July, 2002 August, and 2003 September, 
respectively. We used the Gecko echellette spectrograph fiber fed by CAFE (CAssegrain 
Fiber Environment) from the Cassegrain to Coude focus (Baudrand & Vitry 2000). Spectra 
were centered at 3947 A which was isolated by a UV grism (300 lines mm" with ~ 60 A 
intercepted by the CCD. The dispersion was 0.0136 A pixer^ and the 2.64-pixel FWHM 
of the thorium-argon (Th/Ar) lines corresponded to a spectral resolution of R = 110,000. 
The detector was a back-illuminated EEV CCD (13.5 /xm^ pixels, 200 x 4500 pixels) with 
spectral dispersion along the rows of the device. 

To remove the baseline from each observation, the appropriate mean darks were sub- 
tracted from all the exposures. Flat-fields were then normalized to a mean value of unity 
along each row. All the stellar exposures observed on a given night were combined into a 
mean spectrum to define a single aperture for the extraction of all stellar and comparison 
exposures, including subtraction of residual background between spectral orders (prior to 
fiat- fielding) . This aperture was ultimately used to extract one-dimensional spectra of the 
individual stellar and comparison exposiues and of the mean, normalized flat-field. This 
one-dimensional flat used to obtain the most consistent flat-flelded spectra possible. All 
the data were processed with standard IRAF (Image Reduction and Analysis Facility) rou- 
tines.^ Wavelength calibration was done using the Th/Ar arcs taken before and after each 
spectrum. We required frequent arcs in order to track the CCD drift or any creaking in the 
system throughout the night. Heliocentric and differential radial velocity corrections were 
applied to each stellar spectrum using IRAF's rvcorrect routine. A specimen, fiat-fielded 
spectrum of v And is shown in Figure 1. The Ca H H & K reversals are weak yet visible 
at 3968 and 3933 A. The final spectra were of high S/N reaching 500 per pixel (or 4300 
A~^) in the continuum and 150 (or 1290 A~^) in the H & K cores. 

Spectra with comparable S/N were taken of two stars known not to have close-in giant 
planets, r Ceti and the Sun (sky spectra were taken at dusk). Table 2 of Walker et al. (2003a) 
lists the five CFHT program stars plus r Ceti, including their U magnitudes, exposure times, 
and typical S/N. 



-'^IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the As- 
sociation of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the 
National Science Foundation. 
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2.1.1. Precise Differential Radial Velocities 

Differential radial velocities (ARVs) were estimated with the fxcor routine in IRAF 
(version PC-IRAF V2.f2) which performs a Fourier cross-correlation on dispersion-corrected 
spectra. We used the first spectrum in the series for each star as the template, hence, all 
ARVs are relative to the first spectrum on the first night of the run. Both the template and 
the input spectrum were normalized with a low order polynomial. The correlation used only 
the ~ 20 A region between the H & K lines, that part of the spectrum bounded by (and 
including) the two strong Al I hues (3942 — 3963 A) as shown in Figure 1. 

The ARVs measured during the 2002 July run for the five '51 Peg' stars can be found 
in Walker et al. (2003a). The ARVs from 2003 September are plotted in Figure 2. After 
the orbit is removed, the average aRv is 17 m s~^. The two stars observed at the highest 
S/N have a^v of 7 and 9 m s~^. This precision may be the best achievable for a single 
spectrum with this spectrograph. The excellent ARVs yield current orbital ephemerides and 
hence accurate phases (± 0.02) for each observation. Using the ephemerides of the planets' 
discovery orbits, we tabulated the 2003 September times of sub-planetary position (0 = 0) 
with revised orbital periods. These are listed in Table 2. 

2.2. VLT Observations 

We obtained high-resolution spectra through Visitor mode using the VLT's Ultraviolet 
and Visual Echelle Spectrograph (UVES) mounted on the 8.2-m Kueyen (UT2) over four 
photometric half-nights (2004 April 4 — 7). The standard blue arm setting was used, centered 
on 4370 A, giving a wavelength range of 3750 to 4990 A. We used the CD2 cross-disperser 



Table 2. 2003 September Ephemerides 



Star 




HJD at^ = 


^(HJD)'* 


Revised Port 


SiPorbT 




m 


days 


days 


days 


days 


T Boo 


33 


2452892.864 


0.066 


3.31250 


0.00026 


HD 179949 


19 


2452894.114 


0.062 


3.09246 


0.00031 


HD 209458 


17 


2452893.653 


0.070 


3.52490 


0.00020 


51 Peg 


15 


2452895.868 


0.085 


4.23092 


0.00014 


V And 


9 


2452892.615 


0.092 


4.61750 


0.00052 



'^Uncertainties in the respective measurements 
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grating (660 g mm-^) with a CCD of 2048 x 3000 pixels of 15 /xm^. Image Slicer #2 with 
a sht width of 0.44" resulted in a resolution of R f^i 75,000. 

The data was reduced on-site by the UVES data reduction pipehne that uses the ESO- 
MIDAS software package within the UVES content. The data processing consisted of stan- 
dard procedures: bias subtraction, interorder background correction, cosmic ray hit removal, 
fiat- fielding, and wavelength calibration. The data were wavelength calibrated with Th/Ar 
arcs attached at the beginning and end of each set of 8 to 10 stellar exposures. The exposure 
times, number of exposures and the S/N for each star are listed in Table 3. 

3. A Compcirison of Ca II Emission 

3.1. Extracting the H, K and Al I Lines 

The very strong Ca II H and K photospheric lines suppress the stellar continuum making 
it difficult to normalize each 60-A spectrum consistently. For this reason, a careful analysis 
by which to isolate any modulated Ca II emission was devised. Figure 1 shows a flat-fielded 
spectrum of v And with the normalization levels marked with dashed lines. The normal- 
ization wavelengths were constant for all spectra of a given star. The 7-A spectral ranges, 
centered on the H and K lines, were chosen to isolate the H and K reversals while minimiz- 
ing any apparent continuum differences induced by varying illumination of the CCD. This 
window is, however, wide enough such that a few photospheric absorption features appear 
which could be tested for variability as well. The mean Ca II K cores for the program stars 
are shown in Figures 3 and 4. Also, 7-A and 2-A cuts, centered on the strong photospheric 
Al I fine at 3944 A were used as internal standards since the line has comparable depth and 
S/N to the Ca II lines. 

To normalize each sub-spectrum, the end points were set to 1 and fitted with a straight 
line. The spectra were grouped by date and a nightly mean was computed for each of the 
three fines. The RMS of the Al I residuals for each CFHT target star is less than 0.0005 
of the normafized mean. This is representative of the Al I residuals for all stars in all four 
runs. The Al I fine of the VLT data varies with an average RMS of 0.0006. These values 
demonstrate both the level of stellar photospheric stability as well as the reliability of the 
data reduction and analysis. For all stars observed, both the H & K emission was non- varying 
at the 0.001 level on a given night, likely as a result of intra-night (statistical) noise. 
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3.2. Long-term Variability 

With four CFHT observing runs spanning a baseline of over two years, we can compare 
the long-term variations in the chromospheric levels of the stars. We measure emission 
strength by integrating across the normalized K cores bounded by the Kl features (Montes 
et al. 1994). Each star's average integrated K emission (K) for each observing run is plotted 
in Figure 5 along with its fractional variation relative to the overall average emission. This 
is a good start to tracking the intrinsic stellar activity cycles of these stars. In the case 
of the Sun, we see the decrease from 2001 August to 2003 September as it declines from 
solar maximum. However, we also observed the naked-eye sunspot grouping of 2002 August 
that appeared in our data as a ~ 2% increase in the Ca II emission relative to the other 
years. Since the variability from run-to-run may also be an indication of active regions on 
the disk of the star, we require more frequent monitoring over several more years to firmly 
say anything about the activity cycle of any individual program star. 

3.3. Short-term Vciriability 

When monitoring Ca II emission, intrinsic stellar activity modulated by stellar rotation 
will appear along with any possible activity stimulated by the planet. The orbital periods 
of the planets are well known and uniquely established by the PRV and transit discovery 
methods. The rotation periods of the stars are much harder to determine in part due to stellar 
differential rotation which yields non-unique periods. For our work, it is key to distinguish 
between the rotational and orbital modulation of chromospheric emission. 

To isolate the chromospheric activity within the reversals, we took nightly residuals 
from the average stellar spectrum of all data. Each residual spectrum had a broad, low- 
order curvature removed which was an order of magnitude less than the variations in the H 
and K hnes discussed below. The residuals of the normalized spectra (smoothed by 21 pixels) 
were used to compute the Mean Absolute Deviation (MAD = N'^Tjldatai — mean\ for N 
spectra). As an example, the nightly residuals used to generate the MAD plot for v And are 
displayed in Figure 6 (top). The MAD plot with the corresponding K-core superimposed is 
in Figure 6 (bottom). The identical analysis described above was performed on the Ca II H 
line of all target stars. The complete Ca II H results can be found in Shkolnik (2004). As the 
two resonance lines share the same upper level and connect to the ground state, the same 
activity is seen in both. As expected, the Ca II H emission and activity levels are 2/3 
that of Ca II K (Sanford & Wilson 1939). 

The star with the shortest rotation period in our sample is r Boo. It has the largest 
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vsini {— 14.8 m s ^; Gray 1982) and is believed to be in synchronous rotation with its 
tightly-orbiting planet {Prot = 3.2 ± 0.5 d, Henry et al. 2000; P^rb = 3.31250 ± 0.00026 
d, this work). If r Boo is tidally locked to its planet, the planet-star interaction may be 
minimal (Saar et al. 2003, 2004) due to the fact that there is near zero relative velocity. (See 
Section 4 for theoretical discussion.) The integrated residuals from the mean normalized K 
core are plotted against orbital phase in Figure 7 (left). Looking at the individual observing 
runs in the plot for r Boo, the two nights of observation in 2001 August did not show 
much variation when the star was somewhat less active. The 2002 July data showed an 
increase in activity near the sub-planetary point relative to the other observations in their 
respective runs. However, in 2002 August and 2003 September observations show a relative 
enhancement at (f) ~ 0.4 — 0.5 modulated with a period near Prot and Porb- This chaotic 
activity does not allow us to draw any conclusions about planet-induced heating. 

Similar to r Boo, the K emission of HD 209458 showed night-to-night modulation, but 
with a smaller amplitude during most runs and without any phase coherence, as shown 
in Figure 7 (right). In 2001 August, we caught the system immediately after transit at 
which time we observed a slight enhancement in the Ca II emission relative to all other 
observations. In the 2002 July run, an increase in emission occurred at ~ 0.25 with no 
apparent rise toward = 0. Due to the relatively low S/N of these data and the large 
intra-night deviations, we cannot form any conclusions. 

As seen in Figures 7 and 8, four of the five CFHT stars show significant chromospheric 
variation throughout a single observing run. The standards, r Ceti and the Sun, show no 
such modulation. One consistent result for the 'active' stars (r Boo, HD 179949, HD 209458 
and V And) is their night-to-night modulation of H & K emission. Also, unhke the case 
for HD 73256 (see Section 3.4), the night-to-night variations of these stars do not increase 
or decrease monotonically throughout an observing run, implying the variability cannot be 
explained exclusively with starspots rotating into or out of view. Another mechanism is 
necessary. The night-to-night variations may indicate planet-induced activity or sporadic 
fiaring from hotspots. If coupled to the planet, the localized activity would be travelling on 
the stellar surface faster than the star is rotating as it tracks the planet in its orbit. 

Other than for r Boo, the timescale of activity is short compared to the stellar rotation 
period. Unfortunately, due to the large uncertainties in the rotation periods of the other 
stars, phasing with rotation is uninformative at this stage. However, we do know that v And 
and HD 209458 both have Prot > ^Porb and rotate only < 20° per day. Wolf & Harmanec 
(private communication) recently made UBV photometric observations of HD 179949 at our 
request from the SAAO Sutherland Observatory. When they combined their V observations 
with those from Hipparcos (converted to V) they detected a rotation period of 7.07 d but 
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with an amplitude of only 0.008 mag. Given that the RMS of the V-observations was 0.006 
mag, this periodicity is at the limit of detection. Indirect indications of the rotation rate 
of HD 179949 imply P^ot ~ 9 days and are presented in Shkolnik et al. (2003). (See also 
Saar et al. 2004.) These include a high X-ray luminosity for the star, a very long tidal 
synchronization timescale and a moderate Shk index. While more photometry is needed 
to determine a rotation period conclusively, it is highly unlikely that HD 179949 is tidally 
locked to its planet at 3.092 d. 



3.3.1. Planet-induced Activity on Two Stars 

In Shkolnik et al. (2003) we presented the first evidence of planet-induced heating on 
HD 179949. The effect lasted for over a year and peaked only once per orbit, suggesting a 
magnetic interaction. We fitted a truncated, best-fit sine curve with P = Porb = 3.092 d 
corresponding to the change in projected area of a bright spot on the stellar surface before 
being occulted by the stellar limb. Figure 8 (left) updates the integrated K residuals to 
include the 2003 September data. The spot model is a remarkable fit for the 2001 and 2002 
data peaking at (f) — 0.83 with an amplitude of 0.027. Clearly, the average K emission is 
higher during the latest run (as shown in Figure 5) with a much smaller level of variability. 
It is interesting to note that the 2003 data still peak between <p — 0.80 — 0.95, consistent 
with the previous results. 

The second convincing case of magnetic interaction is between v And and its innermost 
giant planet.^ In Figure 8 (right), the 2002 July, 2002 August and 2003 September runs show 
good agreement in phase-dependent activity with an enhancement at = 0.53. The best-fit 
sine curve has an amplitude of 0.0032. The 2001 August fiuxes are lower than the mean of 
all four observing runs by almost 3% and still display a significant (> 2a) modulation like 
the quiescent epoch of HD 179949. Again, even the low-amplitude modulation has a rise 
and fall with a period consistent with P(yrh and peaks near =0.5. 

For these two cases, the peak of the emission does not directly coincide with the sub- 
planetary point, = 0. For HD 179949, it leads the planet by 60° in phase and for v And, 
the Ca II emission is 169° out of phase with the sub-planetary point. Santos et al. (2003) 
also observed a 90° lag from the sub-binary point in the periodic activity indicated by the 
photometric variations for HD 192263. The phase lead or lag may help identify the nature 
of the interaction. For example, the phase offset of a starspot or group of starspots can be 
a characteristic effect of tidal friction, magnetic drag or reconnection with off-center stellar 



'^v And has three known Jupiter-mass planets at 0.059, 0.829, and 2.53 AU (Butler et al. 1999). 
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magnetic field fines, including a Parker spiral-type scenario (Weber & Davis 1967, Saar et 
al. 2004). In any case, the phasing, amplitude and period of the activity have persisted for 
over a year between observations. For HD 179949, this equals 108 orbits or at least 37 stellar 
rotations and for v And, the time spans 88 orbits or approximately 29 rotations. 

The observations are consistent with a magnetic heating scenario as the chromospheric 
enhancement occurred only once per orbit. We estimated the excess absolute flux released 
in the enhanced chromospheric emission of HD 179949 by calibrating the flux with that of 
the Sun. The flux was the same order of magnitude as a typical solar flare, ~ 10^^ erg s^^ or 
1.5x10^ ergs cm~^ s~^. This implies that flarc-likc activity triggered by the interaction of a 
star with its hot jupiter may be an important energy source in the stellar outer atmosphere. 
This also offers a mechanism for short-term chromospheric activity on the stars with close-in 
Jupiter-mass planets. 



3.3.2. The Non- Varying Program Stars 

Of the 10 program stars wc monitored for H & K variability, five of them showed no 
changes down to the 0.001 level: 51 Peg, HD 46375, HD 75289, HD 76700, HD 83443. There 
are two reasons we offer to explain the relative quiescence of these stars. It was well-known 
through the many years of the Mt. Wilson Shk survey there is a strong correlation between 
rotation rate (or inversely with rotation period) and Ca II emission (Noyes et al. 1984, 
Pasquini et al. 2000). This is a likely contributor though is not obviously clear in our sample 
set as shown in Figure 9 (left) where the inverse of the rotation period is plotted against 
the mean chromospheric emission of the Ca 11 K line (K'). The photospheric contribution 
to the emission was removed from (K) using (K') = (K)(l - -^^) where Rphot = ^hk — 
R'hk and is an empirical function of {B — V) and Shk taken from Hartmann et al. (1984). 
The chromospheric contibutions R'hk are from Wright et al. (2004).^ In Table 1 we list (K), 
(K'), and (MADK) for all the stars. From Figure 9 (right), we deduce that the higher a star's 
chromospheric emission, the more night-to-night activity it displays. Radick et al. (1998) 
show the same effect for a much larger sample. This is akin to shot noise since the fiaring 
or stochastic noise associated with the activity will increase with the activity level. 

Secondly, a recent calculation of the magnetic fields in giant extrasolar planets (Sanchez- 
Lavega 2004) looked at the internal structure and the convective motions of these planets 
in order to calculate the dynamo-generated surface magnetism. Given the same angular 



'''For those few stars that were not in Wright et aL's paper, wc removed tlie photospheric contribution as 
tabulated from stars in Wright et al.'s sample of the same spectral type and log((?). 
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frequency (which is a reasonable approximation for the short-period planets in question) , the 
magnetic dipole moment, and hence the magnetospheric strength, increases with planetary 
mass. This is observed in our own solar system for the magnetized planets where the magnetic 
moment grows proportionally with the mass of the planet (Argc ct al. 1995). Since only lower 
limits exist for most of the hot jupiters, we can only plot MpSim against (MADK) in Figure 10 
where we still see an intriguing correlation. The dashed circles for HD 179949 and v And are 
their (MADK) values (0.0021, 0.0011, respectively) with the orbital modulation removed. 
Of our sample, r Boo has the most massive planet and yet falls well below the correlation. 
As we discuss further in Section 4, if the star and planet are tidally locked, as is thought to 
be the case for r Boo, then there is httle or no free energy left from the orbit and we would 
expect weak, if any, magnetic coupling. 



3.4. Modulation by Rotation 

For most of our target stars, rotation periods are not well enough known to accurately 
phase the Ca II data with P^ot- There is even ambiguity in Prot of the often-observed k} Ceti. 
In late 2003, k} Ceti was monitored continuously for 30.5 days by MOST.^ This best-ever 
lightcurve obtained for Ceti showed a pattern composed of two transiting spots of differing 
periods, 8.9 and 9.3 days providing a direct and unique measurement of differential rotation 
(Rucinski et al. 2004). The real ambiguity in P^-ot from the space observations of k} Ceti 
demonstrates the difficulty in measuring P^ot for most stars. 

Nonetheless, two stars do exhibit clear rotational modulation: Ceti, a star with no 
confirmed planet (Halbwachs et al. 2003) and HD 73256, a star with a 1.85-Mj planet orbiting 
at 0.037 AU (Udry et al. 2003). The Ca 11 emission from Ceti has been monitored by 
Baliunas et al. (1995) through the narrow-band filter of the Mt. Wilson survey from 1967 to 
1991. These data show long-term stability of a period of 9.4 ± 0.1 d, close to the photometric 
rotation period of 9.214 d pubhshed by Messina & Guinan (2002). The rotation period for 
HD 73256 is photometrically determined to be 14 d (Udry et al. 2003), consistent with the 
13.9 days derived from the ^ hk activity index (Donahue 1993). 

We observed periodic Ca II H & K variabihty in the chromosphere of Ceti during 
our 2002 and 2003 CFHT runs from which we determined a rotation period of 9.332 ± 0.035 
d. The results were first pubhshed in Rucinski et al. (2004) where we compared the activity 
seen in Ca II with MOST's lightcurve showing that the chromospheric activity coincided 



^Thc MOST microsatellite (Microvar lability & Oscillations of STars), a Canadian photometric telescope 
recently launched to observe p-mode oscillations on sun-like stars (Walker et al. 2003b) 
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with the low-latitude spot such that the maxima of the two curves agree. HD 73256 was 
observed twice per night for four nights at the VLT. The mean K cores of these two stars arc 
shown at the top of Figure 4 where their similarly strong emission is evident. The high level 
of chromospheric emission points to a young age for both of these stars: 650 — 750 Myr for 
Ceti (Giidel et al. 1997; Dorren & Guinan 1994) and 830 Myr for HD 73256 (Donahue 
1993). The (MADK) given for these two stars in Table 1 are corrected for the rotational 
modulation resulting in activity levels comparable to HD 179949. The integrated residual K 
fluxes for k,^ Ceti and HD 73256 are plotted against relative rotational phase in Figure 11. 
For completeness, we plot the residual K fluxes of HD 73256 as a function of orbital phase 
in Figure 12 where the flare is apparent at 0or6 — 0.03. There is no clear signature of the 
planet in the activity. 

In both Ceti (My = 4.92) and HD 73256 {My = 5.27) show sporadic flaring 

beyond the clear rotational modulation. The periodic best-fit curve for k} Ceti necessitates 
a non-zero eccentricity while a sine curve is sufficient for HD 73256. The largest excursion 
from the rotation curve of Ceti has an energy of 2.8 x 10'^ erg cm^^ s^^ (again, measured by 
comparing with solar absolute flux). We estimate the absolute flux emitted from HD 73256's 
flare at (prot = 0.15 to be 4.9x10^ erg cm~^ s~^ (or > 2.9 x 10^ erg cm~^ if the flare lasted 
for at least the hour for which we observed it.) The modulation of the K emission due to 
rotation is ~ 6% indicating that the emission is dominated by a large hotspot on the stellar 
surface. As we have seen on HD 179949 and v And, planet-induced variations are at the 
level of 1 — 2% suggesting that the reigning hotspot could have diluted any heating caused 
by the hot jupiter. 

The ARVs for k} Ceti are plotted in Figure 13 against the 9.332-d phase determined 
from the K-line residuals where the open symbols correspond to data from 2002 and the 
solid squares to 2003. From the 2002 data we derive aRv =21.8 m s~^ and 23.6 m s~^ when 
combined with 2003. These values are very similar to any —'^4:A m s~^, found over 11 years 
by Cumming et al. (1999). In 2003 the ARVs appear significantly different between the 
two nights, something which seems to be reinforced by the consistency within the pairs of 
ARVs. While a planetary perturbation cannot be ruled out by the 2002 data and other PRV 
studies, the difference in 2003 might be associated with the velocity field of the star itself. 
The increase of velocity with increasing K-line strength at (prot ~ 0.7 is consistent with the 
extreme event in 1988 seen by Walker et al. (1995). However it should be emphasized, the 
possibility remains of a close giant planet around Ceti. For instance, a planet inducing 
a reflex radial velocity variation < 50 m s~^ and tidally synchronized with the star, would 
have MpSini ~ 0.74 Mj and a = 0.084 AU. 
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4. A Physical Scenario 

The enhancements of chromospheric activity on HD 179949 and v And appear only 
once per orbit implying a magnetic, rather than tidal, interaction between the star and its 
hot jupiter. The two stars are F-type stars with higher X-ray luminosity than the solar 
value. The ROSAT catalogue of bright main-sequence stars hsts HD 179949 as having at 
least double the X-ray luminosity (a measurement independent of i) of most other single 
F8 — 9 dwarfs (Hiinsch et al. 1998). Yohkoh solar X-ray observations (e.g. Yokoyama 
& Shibata 2001) have shown that the energy release at the site of magnetic reconnection 
during a solar flare generates a burst of X-ray emitting gas. This hot plasma is funneled 
along the magnetic field lines down to the surface producing 'footprints' through anomalous 
heating of the Sun's chromosphere and transition region. It is this same phenomenon that is 
likely occurring between hot jupiters and their host stars through the reconnection of their 
magnetic fields. Observationally, companion-induced activity is unambiguously observed on 
RS CVn stars, as discussed in Section 1. 

The hot jupiter of v And is located farther out from the host star than that of HD 179949, 
implying that the magnetic interaction between the star and its hot jupiter is diminished. 
This would result in less Ca II enhancement as is shown in our data, r Boo is also an F-type 
star with intense X-ray emission, but its magnetic infiuence is likely reduced by the small 
relative motion in the azimuthal direction between the planet and the stellar magnetosphere 
in an almost final equilibrium state in which the star and the planet are tidally locked to 
each other. While more data are required to truly verify whether the variability of Ca II 
emission is correlated to the apparent position of the hot jupiter in these systems, we explore 
and review a few theoretical aspects of planet-induced heating scenarios in this section. 

The location of the hot jupiter relative to the Alfven radius (the distance from the star 
at which the radial velocity of the wind K-,wind equals the local Alfven velocity Va) plays a 
significant role in transporting energy toward the star against the stellar winds. Since the 
Alfven radius of the Sun is about 10 to 20 R0 at solar minimum and 30 R0 at solar maximum 
(e.g. Lotova et al. 1985), the small distance < 0.1 AU of hot jupiters from their host stars 
suggests that unlike our Jupiter, surrounded by a bow shock, some of these hot jupiters 
are located inside the Alfven radius depending on the magnetic strength of their host stars 
(Zarka et al. 2001, Ip et al. 2004). Therefore the direct magnetic interactions between a hot 
jupiter and its star without a bow shock might resemble the lo- Jupiter interactions (Zarka 
et al. 2001) or the RS CVn binaries (Rubenstein & Schaefer 2000). Most of the theoretical 
models applied to these two cases have focused on the geometry of intertwined magnetic 
fields as well as the the energy transport through Alfven waves and/or induced currents. 

Alfven waves cannot propagate along the stellar field lines toward the star in the region 
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outside the Alfven radius where the group velocity of Alfven waves is always in the positive 
radial direction (e.g. Weber & Davis 1967).^ Other means of inward energy transport require 
their energy flux to be at least larger than the energy flux carried by the stellar winds. 

Ip et al. (2004) estimate the input magnetic power due to the relative motion between 
the synchronized hot jupiter and the stellar magnetosphere to be ~ 10^^ ergs s~^, the same 
order of magnitude of a typical solar flare. A similar amount of power might be obtained 
based on the induced current model (Zarka ct al. 2001) if the radius of the ionosphere of an 
unmagnetizcd hot jupiter in the induced current model can be approximated by the radius 
of the magnctopausc of a magnetized hot jupiter. The observed excess energy flux from an 
unresolved disk of the star is equal to this input magnetic power averaged over the disk of 
the star. That is, the energy flux is roughly equal to 

{Bl/ST,){Vorb - V^^,wind)(r^/i?*)' 

\^ 87r )''\Porb Prot) U/ UJ ' 

where r„ is the radius of the planet's magnetopause, is the magnetic field at the magne- 
topause, is the radius of the star, Rp is the radius of the planet, a is the distance between 
the star and the planet, Vorh and Porh are the orbital velocity and period respectively, Prot is 
the rotation period of the star, V^^^ind is the azimuthal component of the stellar wind veloc- 
ity, and and Bp are the mean magnetic field on the surface of the star and the planet. In 
deriving the above equation, we have assumed that the stellar magnetic field decays as r~^, 
the planetary field decays as r"*^, was determined by equating the stellar and planetary 
fields at the magnetopause (i.e. B^{R^/aY = Bp{Rp/rm)'^)i and the stellar magnetosphere 
inside the Alfven radius is nearly co-rotating with the stellar rotation (i.e. V^^wmd ~ aVt^). 
The radial component of the stellar wind K-,wind is left out from the above estimate as long 
as the condition K,wind i$ Krb ~ V'^.wind is valid. The energy flux in eq(l) should be regarded 
as the maximal input energy from the planet's orbital energy because only some fraction of 
this amount of energy is transferred to the Ca II emissions. For p — 2 (Vrsnak et al. 2002 
for the case of our Sun, and Weber & Davis 1967 for the case of open fields), q — 2, a — 
0.045 AU, R^ = 1.3 R0, Rp = 1.1 Rj, Prot = 9 d, = 200 G, and Bp = 10 G, the energy flux 



^Since the sound speed is much less than the Alfven speed in the extended coronal region, the fast Alfven 
radius is almost the same as the Alfven radius. Therefore we do not distinguish these two radii in this 
paper. However unlike the Alfven modes channeling along field lines, the fast Alfven waves can propagate 
in all directions inside the Alfven radius and therefore their impact on the parent star is attenuated. Here 
we assume that the Alfven mode is the dominant, or at least the comparable, mode compared to the other 
compressible modes when they are excited. 
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given by eq(l) is roughly equal to 10^ erg cm~^ s""^, a value comparable to the differential 
intensity from our data of Ca II K emission from HD 179949. The same amount of energy 
flux can be also achieved for the case of a dipole field for the hot jupitcr {p = 2, q = 3) where 
a = 0.045 AU, = 250 G, and Bp = 10 G. If both the star and the planet have dipole 
fields {p = q = S), very strong fields 5* = 1000 G and Bp = 30 G are required to generate 
the same energy flux. Therefore the tight energy budget constrained by the synchronous Ca 
II emission from HD 179949 strongly suggests that the mean global fields of this F-type star 
are not likely in a dipole field configuration at the location of the planet but have the radial, 
open structure (i.e. p — 2) just like the solar fields have as a result of the outfiowing winds. 

The argument that hot jupiters might have weaker fields than our Jupiter due to slower 
spin rates and weaker convection (Sanchez-Lavega 2004) should be treated with caution 
since in addition to the uncertainty in the interior structure of the metallic hydrogen region 
of a hot jupiter, the response of slow convection to various rotation rates in the dynamo 
process is not well understood. If p = 2, and -Bp < 1 G (Sanchez-Lavega 2004) for HD 
179949 b, then eq(l) indicates that B^ = 300 G is required to generate the energy flux 10^ 
erg cm~^ s~^ and in this case the stellar field dominates over the planetary field even on the 
surface of the hot jupiter (i.e. g = oo in eq(l)). The observational constraints on the strength 
of the stellar magnetic fields such as the field versus Ca II relation (Schrijver et al. 1989, 
1992) along with the radio cyclotron emissions from the hot jupiter's magnetosphere^ should 
help to narrow down the field strength of hot jupiters in the magnetic interaction scenario, 
therefore improving our knowledge of the interior structure and the dynamo processes of 
gaseous planets. 

Now we turn our attention to the variation of Ca II level, the phase coverage of additional 
emission, and the phase lead at the different epochs of our observations. The data for 
HD 179949 seem to suggest that the additional Ca II emission is smaller and the range of 
phase spanned by the emission is larger when the average (or minimal) K emission is larger. 
Presumably this has something to do with the intrinsic stellar activity and therefore the 
stellar field geometry. The observed phase lead may be caused by spiral stellar fields loaded 
with stellar winds. When a star is in its quiet phase and therefore the average Ca II emission 
is very low, no enhancement happens probably because Vwmd ^ for a hot jupiter located 
far out from its parent star and therefore being outside the Alfven radius, as may be the case 
for the 2001 August observations of v And. As the star shifts to its more active phase, the 
average Ca II emission is at the moderate level and V^ind is not far smaller than Va at the 
hot jupiter. At this time, the configuration of the field lines on the surface of the star may 



^Thc characteristic frequency of the cyclotron radiation is c ' ^^^^^ ^ the speed of hght, e and rrie 
are the electron charge and mass respectively. The radiation can reach w 30 — 60 MHz if w 10 — 20 G. 
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be characterized by the open fields from the coronal holes covering a large area of the star, 
as well as the closed fields distributed only near the magnetic equator. Consequently the 
foot-points of the open stellar field lines pointing to the planet arc located at lower latitudes, 
as indicated by the 2001 and 2002 data for HD 179949 fitted by a truncated sine curve 
(Shkolnik et al. 2003; also shown in Figure 8). While a great number of closed fields form 
at low magnetic latitudes during this time, the coronal holes shrink to the magnetic polar 
regions. The bright spot at a high latitude imphed from the observation for HD 179949 in 
September 2003 might be caused by the scenario that the hot jupiter perturbs the open field 
lines emanating from the shrinking coronal holes near the polar regions, leading to a longer 
phase duration of the additional Ca II fiux and perhaps reducing the additional emission due 
to the smaller projection along the line of sight at higher latitudes. Note that at the time 
close to the stellar-maximum activity, the stellar field lines might be occasionally stretched 
out like solar streamers emanating from the low latitudes of the star where closed loops of 
stellar magnetic fields aggregate, possibly giving rise to the planet-induced heating at low 
latitudes of the stellar surface as well. 

The picture that we have sketched thus far assumes that most of energy fiux released 
from the vicinity of the hot jupiter is transported along the field lines by Alfven waves 
and deposited at the foot-points of the magnetic lines. Since the field lines inside the Alfven 
radius are dominated by the poloidal component, detailed calculations for stellar wind models 
are needed to study how the integration of small pitch angles of the field fine can lead to 
the moderate to large phase angles, 60° for HD 179949 and 180° for v And. The phase 
difference is also determined by the Alfven-speed travel time along the stellar field line. If 
Va» Vwind- the Alfven disturbance roughly takes a/V^ ~ a few hours to propagate from 
the hot jupiter to the star with Va ~ lO'^ — 10^ cm s^^ at 0.04 AU. This means that the 
azimuthal angle that the planet has already traveled over the Alfven-speed travel time is 
not small. Since the hot jupiter of v And is located farther out (a = 0.06 AU) from its 
star than HD 179949 b (a = 0.045 AU), the large phase angle of 180° for v And might 
actually represent a phase lag caused by the small inward group velocity of the Alfven waves 
Va — K,wind that takes considerable amount of time to travel along the field fines right after 
the waves are generated from the planet. Alternatively, a large phase difference between the 
location of the heating spot in the chromosphere and the position of the hot jupiter might be 
caused by the entangled or rotationally spiraled magnetic fields connecting directly between 
the hot jupiter and the star (Rubenstein & Schaefer 2000, Saar et al. 2004, Ip et al. 2004). 

Besides considering the energy transported along the field lines and field geometry, 
understanding how the energy is dissipated is important to construct a complete picture for 
the Ca 11 emission in the planet-induced scenario. In our own solar system, a bow shock 
of the solar wind generally hinders this because Alfven waves cannot propagate toward the 
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star from the region outside the Alfven radius. The inward Poynting energy of Alfven 
waves r^}?Bj ^ is roughly conserved if the wind velocity K.wind is much smaller than the 
Alfven speed Va- Here r is the distance away from the star, h is the magnetic disturbance 
of the wave, and p is the mass density of the stellar wind. At first glance, the transition 
region characterized by a steep decrease of the Alfven speed with depth seems to provide a 
possible radial stratification for the inward-propagating Alfven waves to pile up the magnetic 
energy density, leading to the nonlinear dissipation of the growing waves and therefore the 
heating on the top of the chromosphere. However, the extremely narrow transition region 
corresponding to sharp gradient in density and Alfven speed should act as a wave barrier 
to reflect the Alfven waves. In this case, the wavelength is comparable to the size of the 
magnctopausc of the hot jupitcr (Wright 1987) unless the high-frequency modes are largely 
excited at the reconnection site. The energy carried by the planet-induced Alfven waves 
along the open fields might be finally transmitted to coronal loops and therefore might be 
dissipated via resonance absorption, a damping mechanism however in contradiction to the 
observations of the solar corona and the coronae of cool stars (Schrijver & Aschwanden 2002; 
Demouhn et al. 2003). The heating due to accelerated particles by Alfven waves (Crary 1997) 
is not important either because unlike the interplanetary space, the dense corona eliminates 
the kinetic effect of plasmas inside 0.04 AU.^ Despite the difficulties, the energy deposit 
on the surface of the star from the planet-induced Alfven disturbances may be achieved by 
interacting non-lincarly with stellar winds. The stellar wind consists of charges particles, 
stellar fields, and probably Alfven waves. Non-linear interactions between planet-induced 
incoming and stellar intrinsic outgoing waves is one route of heating the surface of the star. 

The scenario of magnetic interaction implies the orbital decay of hot jupiters since the 
ultimate source of energy comes from the orbital energy. Theoretically the orbital decay 
of hot jupiters on the timescale of a few billion years can result from the tidal dissipation 
in the host star driven by the hot j up iter so long as the tidal dissipation in the solar-type 
stars is efficient (Rasio et al. 1996, Witte & Savonije 2002, Patzold & Rauer 2002, Jiang et 
al. 2003). In the magnetic interaction scenario for Ca II emission, the time scale of orbital 
decay is roughly equal to the ratio of the orbital energy of the hot jupiter (~10^^ ergs) to 
10^^ erg/s. This gives a timescale as short as several billion years, imposing a non-negligible 
constraint on modeling the orbital evolution of hot jupiters. 



^Besides the particle acceleration by Alfven waves, a jet of accelerated particles with high speeds from the 
reconnection site at the magnetopause of a magnetized planet might be able to heat the stellar corona, but 
it is difficult for the process to produce a large phase difference between the heating spot and the location 
of the planet unless the coronal fields can be entangled in a more large-scale manner like the models for the 
RS CVn stars. 
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In the case of r Boo, the F-type star might be almost tidally locked to its hot jupiter. 

Therefore there is less free energy available from the planet's orbit. According to eq(l), the 
azimuthal relative motion between the stellar and the planetary magnetosphcres is smaller 
than that for the HD 179949 system by a factor of 0.2 if 3.2 d is indeed the rotation period 
of r Boo. The mean chromospheric emission (K') from r Boo is weaker than that from HD 
179949 roughly by a factor 0.9. Assuming that the mean stellar field of r Boo is smaller 
than that of HD 179949 by 0.9^ ((K') oc B°-^), one can estimate that the input energy to r 
Boo is less than that to HD 179949 roughly by a factor of 0.15. If the corresponding Ca II 
emissions are dimmed by this same factor, this effect should have been detected. However 
unhke the other two F-type stars HD 179949 and v And, the variability of Ca II emissions 
from r Boo did not show any consistent phase relation with the planet's orbit. Note that 
the radial movement of spiral stellar fields with the stellar winds might be important in 
providing energy in this case because Kswind ~ Vorb at 0.04 AU in the solar system. However 
in the case of small pitch angles of the field lines inside the Alfvcn radius, the Alfvcn modes 
driven solely by the radial impact of the stellar winds might not be as efficient as the stellar 
fields sweeps across the hot jupiter due to the slow relative azimuthal motion. 

The transiting system HD 209458 did not show synchronized Ca II enhancement, per- 
haps because G stars have K-,wind > Va at the distance of the hot jupiter due to weaker 
stellar fields. The non-synchronized Ca II enhancement of r Boo and HD 209458 (as well as 
the fiaring on HD 73256 and k} Ceti, discussed in Section 3.4) might still be attributed to 
direct interactions with their planets but through more chaotic means due to the variability 
of the stellar winds at the locations of these hot jupiters relative to the nearby Alfven radius. 
Monitoring these stars continuously through several orbital and rotational periods should 
pin down the cause — intrinsic due to fast rotation or induced by a hot jupiter — of the 
strong night-to-night variabilities detected in these systems. 



5. Summary and Future 

Of the sample of stars observed from the CFHT, those with planets (with the exception 
of 51 Peg) show significant night-to-night variations in their Ca II H and K reversals, r Ceti 
and the Sun, which have no close planets, remained very steady throughout each of the four 
observing runs. HD 179949 and v And exhibited repeated orbital phase-dependent activity 
with enhanced emission leading the sub-planetary point by 0.17 and 0.47 in orbital phase, 
respectively. Both systems arc consistent with a magnetic heating scenario and may be 
the first ghmpse at the magnetosphcres of extrasolar planets. The phase-lead or lag of the 
peak emission relative to the sub-planetary longitude can provide information on the field 
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geometries and the nature of the effect such as tidal friction, magnetic drag or reconnection 
with off-center magnetic fields, including a Parker-spiral type scenario. 

r Boo and HD 209458 also exhibited night-to-night variations that could not exclusively 
be due to stellar rotation. If r Boo is indeed tidally locked by its hot jupiter, there is no 
orbital energy available to generate the Alfven modes efficiently as the stellar fields sweep 
across the planet due to the slow relative azimuthal motion. We measured the excess absolute 
flux released in the enhanced chromospheric emission of HD 179949 to be the same order of 
magnitude as a typical solar flare, ~ 10^^ erg s^^ or 1.5x10^ ergs s^^ cm^^. This implies 
that flare-like activity triggered by the interaction of a star with its hot jupiter may be 
an important energy source in the stellar outer atmosphere. This offers a mechanism for 
short-term chromospheric activity. 

The H & K emission of k,^ Ceti, an active star with no confirmed planet, was clearly 
modulated by the star's 9. 3-d rotation. Similarly, HD 73256 displayed rotational modulation 
with its 14-d period. In these two cases, the chromospheric emission increases by ^ 6% 
(relative to the normalization level at 1/3 of the continuum). Any planet- induced heating 
at the level of 1 — 2% could have been diluted by the dominating hotspot on the stellar 
surface. Neither the ARVs nor the Ca II periodicity exclude the possibility of a sub-stellar 
companion in a tight orbit around Ceti. 

Apart from the cyclical component for four of the stars, short-term chromospheric ac- 
tivity appears weakly dependent on the mean K-line reversal intensities for the sample of 13 
stars. Also, a suggestive correlation exists with MpSini and thus with the hot Jupiter's mag- 
netic field strength. Because of their small separation (< 0.1 AU), many of the hot jiipiters 
lie within the Alfven radius of their host stars which allows a direct magnetic interaction 
with the stellar surface. 

Additional Ca II observations are crucial to confirm the stability of the magnetic inter- 
action as well as to establish better phase coverage. Observations on timescales of a few years 
will begin to characterize the long-term activity of our program stars and allow us to see 
correlations between intrinsic Ca II emission and night-to-night activity more clearly. This 
work opens up the possibility of characterizing planet-star interactions with implications for 
extrasolar planet magnetic fields and the energy contribution to stellar atmospheres. 

A next step in understanding planet-star interactions is to map the activity as a function 
of stellar atmospheric height. Above the chromosphere lies the thin transition region (TR), 
where the temperature increases steeply as density and pressure drop, and the corona, which 
can extend out to several stellar radii. And since the magnetic field drops off as (where 
2 < p < 3), these layers facihtate a stronger interaction with the planet. Their FUV and 
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X-ray emissions will be extremely important diagnostics. One indication that the heating 
is from the outside in is if the increase in emission occurs slightly earlier in phase than in 
Ca 11. Moreover, the relative strengths of the different emission lines will tell us where most 
of the energy is dissipated. The energy sum will point out if there are any discrepancies 
with the theorized energy budget. Orbital phase-dependent variability at these heights will 
constrain further the nature, form and strength of the interaction as well as specify non- 
thermal radiative processes in these hot layers of gas. 
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Council of Canada (D.A.B.) is gratefully acknowledged. We are indebted to the CFHT 
staff for their care in setting up the CAFE fiber feed and the Gecko spectrograph, as well 
as to the staff at ESO's VLT for their telescope and instrument support and the real-time 
data-reduction pipeline. Also, we appreciate the helpful comments and suggestions from the 
referee, Steve Saar. 
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Table 3. The VLT Program Stars 



Star 


U 


B 


Exp. time 














s 




cont 


core 


HD 46375 


9.33 


8.7 


300 


8 


510 


150 


HD 73256 




8.86 


300 


9x2 


520 


150 


HD 75289 


7.04 


6.94 


120 


10 


790 


230 


HD 76700 


9.17 


8.76 


300 


8 


540 


160 


HD 83443 


9.54 


9.03 


300 


8 


420 


120 



'^Number of spectra taken per night. We observed HD 73256 
twice per night for a total of 18 exposures. 

''Nightly average per 0.015 A per pixel in the continuum near 
3950 A and in the Ca II K core. 
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Fig. 1. — A single spectrum of v And. The arrow defines the 3942 — 3963 A region used to 
measure the differential radial velocities. The dashed hues show the levels of normahzation 
use to isolate the H and K cores. 
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Fig. 2. — The differential radial velocities for the five '51 Peg' stars measured in 2003 
September are plotted as a function of relative orbital phase, and the unvarying star, r Ceti, 
as a function of time. Each sine curve has the planetary orbital period from Table 1 and 
pubhshed amphtude for the star and has been shifted in phase to give the best fit to the 
ARVs. The measurement error of the individual points is < 5 m s~^. The a^v of the ARV 
are listed in Table 2. Note that r Ceti's ARVs are plotted on a much smaller scale than the 
others and have a a^v of 7 m s~^, making the average a^v for all the data 17 m s~^. 



-32- 




Fig. 3. — The mean normalized Ca II K cores for the CFHT program stars. 
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Fig. 4. — The mean normalized Ca II K cores for the VLT program stars plus k} Ceti, which 
we observed at the CFHT. 
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Fig. 5. — Top: the mean integrated Ca II K emission from each of the four runs for the 
CFHT program stars, including r Ceti and the Sun. Units are in equivalent angstroms 
relative to the normalization points. Bottom: the fractional variation for each star relative 
to its overall mean K emission. 



- 35 - 




0.005 

_o 

I 0.004 
Q 

I 0.003 
o 

c 0.002 



0.001 




3932 



3933 3934 
Wavelength (A) 



3935 



Fig. 6. — Top: The residuals (smoothed by 21 pixels) from the normalized mean spectrum 
of the Ca II K core of v And. Bottom: The mean absolute deviation (MAD) of the Ca II K 
core of V And. The units are intensity as a fraction of the normalization level at 1/3 of the 
continuum. Overlaid (dashed line) is the mean spectrum (scaled down) indicating that the 
activity on v And is confined to the K reversal. 
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Fig. 7. — Integrated flux of the K-line residuals from a normalized mean spectrum of r Boo 
and HD 209458 as a function of orbital phase. The symbols distinguish data from different 
observing runs: circles — 2001 August, squares — 2002 July, triangles — 2002 August, 
diamonds — 2003 September. Units of the integrated flux are in equivalent Angstroms 
relative to the normalization level which is approximately 1/3 of the stellar continuum. The 
error bars in residual flux are ± 1 o" as measured from the intra-night variations. The size 
of the phase error is within the size of the points. 




Fig. 8. — Integrated flux of the K-line residuals from a normalized mean spectrum of 
HD 179949 and v And as a function of orbital phase. The solid lines are best-flt spot 
models. For HD 179949, we used a truncated sine curve fltted to the 2001 and 2002 data 
and for v And, we fltted the curve to the 2002 and 2003 data. 
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Fig. 9. — Left: 1/Prot in days plotted against (K'), the mean chromospheric emission in the 
K hne, for all stars less HD 76700 for which no Pro* is published. Right: (K') as a function 

of the mean MAD K values per run for all the stars. The dashed circles arc (MADK) for 
HD 179949 and v And with the orbital (geometric) modulation removed. Units for (K') 
and (MADK) are in equivalent angstroms relative to the normalization level. All values are 
listed in Table 1. 
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Fig. 10. — The minimum planetary mass (in Jupiter masses) plotted against the average 
MAD of the K-hne per observing run. The dashed circles are (MADK) for HD 179949 and 
V And with the orbital (geometric) modulation removed. All data are listed in Table 1. 
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Fig. 11. — Integrated flux of the K-line residuals from a normalized mean spectrum of 
HD 73256 and Ceti as a function of rotational phase. The sohd hues are the best-fitting 
periodical curves. In the Ceti plot, the open circles are 2002 data and the closed circles 
are 2003 data. 




Fig. 12. — Integrated flux of the K-line residuals from a normalized mean spectrum of 
HD 73256 as a function of orbital phase. 
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Fig. 13. — Relative radial velocities for k} Ceti. Open squares are the 2002 data and the 
filled in squares are the 2003 data. 



